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We perform molecular dynamics simulation on a glass-forming liquid binary mixture with the soft- 
core potential in three dimensions. We investigate crossover of the configuration changes caused by 
stringlike jump motions. With lowering the temperature T, the motions of the particles composing 
strings become larger in sizes and displacements, while those of the particles surrounding strings 
become smaller. Then, the contribution of the latter to time-correlation functions tends to be long- 
lived as T is lowered. As a result, the relaxation time and the viscosity T) grow more steeply than 
the inverse diffusion constant at low T, leading to breakdown of the Stokes-Einstein relation. At 
low T, the diffusion occurs as activation processes and may well be described by short-time analysis 
of rare jump motions with broken bonds and large displacements. Some characteristic features of 
the van Hove self-correlation function arise from escape jumps over high potential barriers. 

PACS numbers: 64.70.Q-,61.20.Lc,66.30.hh,66.30.-h 



I. INTRODUCTION 

Recently, much attention has been paid to the slow dy- 
namics in glass- forming liquids [1 , where the relaxation 
time Tq becomes extremely long and the viscosity 77 grow 
dramatically with lowering the temperature T . In par- 
ticular, it is well known experimentally that the curve 
of log 77 vs 1/r (the Angell plot 2 ) becomes steeper with 
increasing 1 /T deviating from the Arrhenius behavior for 
fragile glass-forming Uquids, while the Arrhenius behav- 
ior has been observed for strong glass formers like silica. 
However, it is not yet well understood how this dynami- 
cal crossover takes place microscopically for fragile glass 
formers. Hence, in this paper, we aim to investigate this 
aspect using molecular dynamics simulation. 

In fragile glass-formers, stringlike motions of mobile 
particles [3H5] play a major role in the structural relax- 
ation. The distribution of string lengths [4] was shown 
to widen with lowering T. However, the behaviors and 
the roles of stringlike motions have not yet been fully 
disclosed in the structural relaxation and the plastic de- 
formations. In this paper, we analyze the configura- 
tion changes with the bond breakage theory [SI [7] , which 
was originally used to detect the dynamic heterogeneity 
O in] and has recently been generalized as a statistical- 
mechanical theory of irreversible particle rearrangements 
[lOj . We then confirm that the configuration changes 
mostly occur as stringlike motions of mobile particles at 
low T. We shall furthermore see that the particles com- 
posing strings and those surrounding them behave dif- 
ferently depending on T in a binary mixture at a fixed 
high density n in three dimensions (3D). With lowering 
T, the former move over longer distances with multiple 
broken bonds, while the latter over shorter distances with 
a single broken bond. The number of the latter is several 
times larger than that of the former because of the large 
coordination numbers in 3D (^ 10). Hence, the crossover 
of the latter motions with lowering T should greatly influ- 
ence the relaxation behavior of time-correlation functions 



and the relaxation time at low T. 

In liquid, the Stokes-Einstein relation Drja/T =const. 
between the diffusion constant D and the viscosity ry has 
been successfully applied even for a microscopic test par- 
ticle diameter a. However, this relation is systematically 
violated in fragile supercooled liquids. [TTJ[T2]. Sillescu 
et aZ.|ll| observed the power law behavior D cx rj^'^ with 
1/ = 0.75 at low T. The origin of this violation has of- 
ten been ascribed to the dynamic heterogeneity without 
clear-cut calculations. A number of simulations have not 
clearly explained this violation [SJ [13 [T5in!5] . We shall 
see that the above-mentioned crossover of the particle 
motions around strings is the main origin of the viola- 
tion for our fragile glass-forming system. 

For sufficiently low T in the glass-forming condition, 
the van Hove self-correlation function Gs{r,t) tends to 
have a minimum at a particle length rm and slowly grows 
in the outer region r > rm forming secondary peaks [T51 - 
[T8| . Here, r represents the particle displacement length 
in a time interval of t. We shall see that its minimum 
value at r = Tin becomes very small at low T, indicating 
fast passage of the particles across potential barriers at 
r ^ r^. These behaviors reflect the fact that the diffu- 
sion arises from escape jumps from well-defined transient 
cages with long life times. Sastry et al. called this the 
landscape-dominated regime |14l . which emerges on ap- 
proaching the glass transition. We shall also see that 
the contribution to the mean square displacement M{t) 
from the particles with large displacements behaves as 
QDt soon after the ballistic time region, where the ther- 
mal vibrational contributions from the other particles in 
the interior r < r,n are removed. 

The organization of this paper is as follows. In Sec. II, 
our simulation method will be explained. In Sec. HI, the 
bond-breakage theory in our previous work [SJ [71 [TU] will 
be used to examine the particle motions around strings 
for various T. In Sec.IV, we will examine the long- 
distance diffusion, which will turn out to be governed 
by the activated dynamics at low T. We will also exam- 
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ine the van Hove self-correlation function Gg (f", t) in the 
landscape-dominated regime. 

II. NUMERICAL METHOD 

In this paper we show results of molecular dynamics 
simulation of binary mixtures composed of two species, 1 
and 2, in 3D at low T . We impose the periodic boundary 
condition without applying shear flow. The composition 
is c = N2/{Ni + N2) ~ 0.5 and the total particle number 
is = iVi + A^2 = 10**. The two species have different 
diameters cti and (72 with cr2/o'i = 1-2. The particles 
interact via the soft-core potential, 

^^a/3W = e(^)''-C„^ (r<r,,t), (2.1) 

where a and /? represent the particle species {a,P = 1,2), 
r is the particle distance, and e is the characteristic in- 
teraction energy. The interaction lengths are defined by 

ac.p^{ao,+ap)/2. (2.2) 

The potential vanishes for r > rcut = 3cri . The constants 
Cai3 ensure the continuity of the potential at r = rcut- 
The masses of the two species satisfy TO2/mi = (i72/o'i)'^. 
The average density is given by n = N/V = OSa^^ as 
in some previous papers [TJ [TOl [HI HH] , where V is the 
system volume. The system length is L = V^^^ = 23.2cti. 
Space and time are measured in units of Ci and 

To = fJiVmi/e. (2.3) 

The temperature T is measured in units of e/ks- 

We started from a liquid state at a high temperature, 
quenched the system to a final low temperature, and 
waited for a long time of order 10^. We imposed a Nose- 
Hoover thermostat in these steps. However, after this 
preparation of initial states, we removed the thermostat 
and integrated the Newton equations under the periodic 
boundary condition in the time range t > 0. Thus, the 
particle numbers, the total volume, and the total energy 
are conserved in our simulation (NVE ensemble). At 
the lowest temperature T = 0.24, the simulation time 
was 5 X 10^ for each run and data were taken in the time 
range 5 x 10"* < i < 5 x 10^, during which we did not 
detect any appreciable aging effects in various quantities 
such as the average potential energy and Fs{q,t) defined 
below. 

We are interested in supercooled states, where the 
structural relaxation time Tq is very long. In terms of 
the self part of the density time-correlation function, 

Fsiq,t) = ^^/exp[zq- Ar,(to,to + i)]\, (2.4) 
j 

Tq, is usually defined at q — 2tt(j^'^ by 

Fs{q,T^) = e-\ (2.5) 



Here, q is the wave vector with q = \q\ and Arj(to, ii) = 
rj{ti) — rj{to) is the displacement vector of particle j. 
In this paper, (• • •) denotes taking the average over the 
initial time to of the time interval and over several sim- 
ulation runs. This is needed for averaging quantities in- 
volving two widely separated times, for which accurate 
results do not follow only by the average over all the par- 
ticles for our system size. 

It is worth noting that the stress time-correlation func- 
tion behaves similarly to Fs{q,t), both considerably de- 
creasing in the early stage due to the thermal motions 
(the so-called f3 relaxation) [5J [12] . Since the time inte- 
gral of the stress time-correlation function is equal to the 
shear viscosity t] in the linear response regime, it is nat- 
ural to expect weak T dependence of the ratio ry/T^. For 
the same soft-core potential (2.1), the following relation 
was previously obtained [15]: 

ry= (27r)-Vi-3A:sTr„, (2.6) 

in the range 1 < Tq. < 10"* at n — O.Scrf^. 

III. BOND BREAKAGE 
A. Stringlike motions and their crossover 

We analyze the configuration changes in the bond 
breakage scheme [S1[71[TU]. They occur as jump motions 
causing irreversible particle rearrangements. 

First, we introduce the concept of bond breakage. At 
an initial time t = to, two particles i € a and j G /3 
(a, /? = 1, 2) are treated to be bonded if 

nj{to) < AiCFap- (3.1) 

Hereafter, (i) — \ri{t) — rj(i)| is the distance between 
these particles at time t. At a later time < = <i > to, this 
bond is treated to be broken if 

r,j{ti)> A2<T^p. (3.2) 

Here, AiUap are slightly larger than the first peak dis- 
tances of the pair correlation functions gapif) and A2 is 
somewhat larger than Ai. In this paper, we set Ai — 1.3 
and A2 = 1.7. At high densities in 3D, the bond num- 
ber around each particle is of order 10 because it is on 
the order of the coordination number. Let Ni,{to) be the 
number of the bonds defined at time to and Nb{to, to + 1) 
be the number of the unbroken bonds at later time to + t. 
Then, the fraction of the unbroken or surviving bonds 
after time t is given by 

F,{t) = {mito,to + t)/N,{to)). (3.3) 

The bond breakage time ti, is determined by [S] [7] 

Fbin)=e-\ (3.4) 

Nonlinear rheology in supercooled states [T] [^U] is gov- 
erned by the bond breakage and is characterized by Tf,. 
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(a)T =0.24, Bi>0 (b)T =0.24, Bi>0 
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FIG. 1: (Color online) (a) B particles (those with Bi = 
Bi{to,ti) > 0) composed of their strings and their neighbors, 
(b) BL particles (those with Bf{to,ti) > 0) composed of 
strings, which have large Bi > 3 due to large displacements. 
These are snapshots at t — ti — to = 10'' = 10~^Ta or at 
(p^{t) = 0.1 for T = 0.24 (see Table 1). The particle col- 
ors represent Bi according to the color bar below. Arrows 
represents Ari{to,ti). 



(a)T=0.28. B^>0 (b)T=0.32, B>0 
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FIG. 2: (Color online) BL particles (those with Bf{to,ti) > 
0) for (a) T = 0.28 and t = ti-to = 450 and for (b) T = 0.32 
and t = ti — to — 50. The corresponding snapshot for T = 
0.24 is given in Fig. lb. Here, t — ti — to is determined from 
(j)^ {t) — 0.1. The particle colors represent Bi according to 
the color bar. The strings become ill-defined with increasing 
T. Arrows represents Ari{to,ti). 



Furthermore, we introduce the broken bond number 
for each particle i in terms of Ai and A2 as [TD| 

Bi{ta,ti) = ^e{Aiaap -r,^{to)) 
j 

x0{nj{ti) ~ A2acf}), (3.5) 

where 0{u) is the step function being zero for u < 
and 1 for m > 0. This number is a nonnegative integer, 
tending to zero as ti — to from Ai < A2 and increas- 
ing to 1,2, •• • upon bond breakage. The particles with 
Bi{to, ti) > may be called B particles, while those with 
l3i{to, ti) = may be called non-B particles. The B par- 
ticles are surrounded by different particle configurations 
at the initial and final times t = t^, and ti. 

The displacements of the B particles may not be large 
if they are neighbors of those undergoing stringlike mo- 
tions. Thus it is convenient to group the B particles 
further into those with large displacement Ar^ > and 
those with small displacement Ar^ < ^,„, where is a 
minimum jump length and 

Ar, = Ar,(to,ti) = |r,(ii) - r,(io)| (3.6) 

represents the displacement in time interval [ioi^i]- 
These particles are called BL and BS particles, respec- 
tively. Their broken bond numbers are written as 

S>(to,ti) =S,(to,ii)e(Ar,(io,ti)-^™), (3.7) 
S<(to,ti) = B,{to,h)0{l-m - An{to,h)). (3.8) 

In this paper, we set im = 0.8. A small change of its 
value does not essentially change our conclusions. See 
SubsecIVD and Fig. 12 on this point. 

In Fig. la, we display only the B particles (about 600) 
a.tt = ti-to = 10" ^ 10-Va for T = 0.24, where the 



particle colors represent Bi — Bi{to,ti). We can see that 
many particles with Bi = 1 surround those with Bi > 2. 
Here, the latter particles have undergone stringiike mo- 
tions with large displacements (> 1). Thus, in Fig. lb, we 
depict only the BL particles with Ar^ > 0.8. Here, the 
average of Ar^ over the BL is 1.29, while that over the 
BS particles is 0.152 (see Fig. 4 below). We also notice 
that a majority of the BL particles (74% here) belong 
to the first (small) species (see Table 1 below). Thus, for 
this time elapse of lO", the BL and BS particles mostly 
form strings and their neighbors, respectively. 

At T = 0.24, a large fraction of the BS particles sur- 
rounding strings transiently moved over relatively large 
distances (> 0.5) at the string formation but eventually 
returned to their original positions (not shown in Fig. lb). 
Though rather rare, we also observed that some particles 
largely moved without string formation and returned to 
their original positions after some time (say, 400) . These 
reversible jumps were observed by Vollmayr-Lee |21| . 

However, with increasing T, the difference between the 
BL and BS particles becomes less distinct. In Fig. 2, we 
elucidate how the strings are changed at T = 0.28 and 
0.32, where the number fraction of BL particles (p^ is 0.1 
(see Eq.(3.10) below for its definition). As compared to 
the snapshot at T — 0.24 in Fig. lb, they are shorter and 
more expanded with smaller Bi, so they should rather 
be treated as clusters. This result is consistent with the 
previous calculation of the distributions of string lengths 
and cluster sizes for various T by Donati et al'^ [S] . With 
increasing T, we also notice that the composition of the 
large particles increases among the BL particles, leading 
to larger volume fractions of the BL particles even at 
fixed (j)> = 0.1 (see Table 1 below). 
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FIG. 3: (Color online) Fractions of the B and BL particles, 
4>b{t) and (j)^{t), respectively, for T = 0.24. They are shown 
for t < 800 (left) and on long times (right). Also displayed 
is the fraction of surviving bonds 1 — Ft{t) nearly coinciding 
with (f)^ {t) on long timescales (right). 
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FIG. 4: (Color onhne) (a) Averages {B>) and {B<) in 
Eqs.(3.13) and (3.14) vs T. (b) Averages (Arsz,) and (Arss) 
in Eqs.(3.15) and (3.16) vs T. Here, 0^(t) =0.1, from which 
t = t\ — to is determined for each T as in Table 1. 



B. Differences between BL and BS particles 

Let us consider the numbers of the B and BL particles 
in the system. The average number fractions of these 
particles are given by 

'^"W = l^T.^OWoM))). (3.9) 



TABLE L Data for five temperatures, where t, (t>b{t), c^(i), 
and c^{t) are those at (f)^ {t) — 0.1. See Fig. 4 also. Here, 
(j>b{t) = (t) + 0b~(i) is the fraction of the B particles, while 
(t) and c^{t) are the compositions of the larger species 
among the BL and BS particles, respectively. 



T 




n 


t 








0.24 


109000 


270000 


10000 


0.38 


0.247 


0.547 


0.28 


1170 


103000 


450 


0.46 


0.309 


0.496 


0.32 


41.0 


749 


50 


0.59 


0.329 


0.492 


0.36 


11.4 


267 


20 


0.61 


0.339 


0.481 


0.40 


5.91 


141 


10 


0.59 


0.330 


0.463 



-^^^W = ^T.((^iB>{to,h))), (3.10) 

j 

where t — ii — to- The (j)b{t) is the fraction of the particles 
with broken bonds, while (f>^{t) is that of the particles 
with broken bonds and large displacements. We may call 
(f>^{t) the string fraction also. In our previous paper [TU]. 
we have introduced a bond-preserving time Tbp by 

1 - 0b(rbp) = e-i, (3.11) 

which is the timescale of each particle to have a broken 
bond. For T = 0.24, we have u ^ 2.7 x 10^ = 2.7r„ and 
Thp = 3.2 X 10'* — 0.32tq. Note that Tf, is longer than Tbp 
by one order of magnitude due to the large coordination 
number in 3D. 

In Fig. 3, (j)b{t) and (f>^{t) are plotted on short and long 
timescales. Salient features are as follows, (i) First, on a 
microscopic timescale (t > 4), 0b (t) quickly approaches 
a small number~ 0.03 due to the thermal motions. (In 
2D, we have obtained an algebraic growth, <j)b{t) ^ t^'^, 
for small t, however [1^.) (ii) Second, in the early stage 
(t < 600) , (j)^ {t) grows linearly as 

cl,>{t)^r-,h. (3.12) 

The coefficient r^^^ is the average frequency of rare jump 
motions per particle in the early stage t <^ Tgt (where 
(t>bii) ^ !)• Therefore, it is analogous to the nucleation 
rate in metastable systems In our case, Tgt is of 

order Tf, in Eq.(3.4). See Fig.5d for the T dependence of 
Tst- (iii) Third, in the whole time range in Fig. 3, 0&(t) 
is much larger than (l)^{t). For t 1, this is because 
several broken bonds are produced around a string as in 
Fig. la. In Fig. 3b, the ratio 4)b{t) / (j)^ {t) is of order 5 for 
t > 10^. In addition, (ty^{t) is close to the broken bond 
fraction 1 — Fb{t), where Fb{t) has appeared in Eq.(3.4). 
They are very close for T < 0.3 in the present case, (iv) 
Fourth, at long times t > 1000, we have 0g'(t) - i" ™. 

In Fig. 4, we show how BL and BS paricles behave dif- 
ferently. Displayed in the left are the averages of Bi{to, ti) 
among the BL and BS particles. They are written as 

(^^) = ^E(^^^(^0'*i))' (3-13) 
<^^) = :^E(^^^(^0'*i))' (3.14) 

where = (pb — 4>^ is the average fraction of the BS 
particles. The average broken-bond number is 3.4 for the 
BL particles and is 1.4 for the BS particles at (j)^ {t) = 
0.1 for T — 0.24. Displayed in the right are the averages 
of the displacements Ar^ in Eq.(3.6) among the BL and 
BS particles, 

(Ar)BL = -^E(^^«^(^»^))' (3-15) 
(Ar)B5 = Ti^E(^^^^(^^^))' (3.16) 
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so their contribution to Fs {q, t) becomes long-lived. 

To examine the second feature in more detail, we define 
the conditional self-time-correlation function, 

F^iq,t) = ^ 5^(0(6,) exp[iq . (Ar,]), (3.17) 



where Bj and Arj are abbreviations of Bj(tQ,ti) and 
Arj{tQ,ti), respectively, and only the B particles are 
picked up in the summation. In this function, the con- 
tribution from the BL particles is less than 1% at any T 
(so Bj may be replaced by B^ in the right hand side of 
Eq.(3.17)). In (a) and (b) of Fig.5, Fsiq,t) and F^{q,t) 
are compared at T = 0.24 and 0.28. In (c), the T depen- 
dence of Fg{q,Ta) is shown at t = Tq. It is about 30% 
of Fs{q,Ta) — eT^ for relatively high T, but it increases 
up to 80% for T = 0.24. This dramatic change serves to 
slow down the relaxation of Fg{q,t) at low T. In (d), we 
display Tq, Tbp, Tst, and Tf,, where Ta indeed grows most 
steeply with lowering T. Due to this crossover, we have 
n/Ta - 20 for Ta < 10-* but n/Ta - 2.7 for t„ - 10^ 
from Table 1. 



IV. DIFFUSION 



FIG. 5: (Color online) Top: F4q,t) in Eq.(2.4) and F^{q,t) 
in Eq.(3.17) for (a) T = 0.24 and (b) T = 0.28. Bottom: 
Fs{q,Ta) vs T in (c) and three characteristic times Ta in 
Eq.(2.5), Tbp in Eq.(3.11), r^t in Eq.(3.12), and n in Eq.(3.4) 
in (d). Here, Ta grows more steeply than the others at low T. 



where Ar^, Bf" , and Bf are abbreviations of Ari{tQ,ti), 
B^{to,ti), and Bf{to,ti), respectively. Remarkably, 
( Ar) B L and ( Ar) bs are distinctly separated and in 
Eqs.(3.7) and (3.8) has bee chosen to be 0.8 between 
them. Table 1 presents more data, where (t) and (t) 
are the compositions of the larger species among the BL 
and BS particles, respectively. The former is the com- 
position within strings and is only 0.25 at T = 0.24. 

Thus the motions of the BL and BS particles (those 
composing strings and those surrounding strings) are in- 
creasingly different with lowering T. This should be a 
characteristic feature of fragile glass-forming liquids [21 
[T21 [H] , where the role of the thermal fluctuations in the 
configuration changes crosses over with varying T. 



C. Steep growth of Ta at low temperatures 

The physical meaning of Tq, is not straightforward. 
First, Fs{q,t) in Eq.(2.4) relaxes to a plateau value 
rapidly due to the thermal motions, which is about 0.8 for 
T = 0.24 in Fig.Sa and about 0.7 for T = 0.28 in Fig.5b. 
Second, with lowering T, the motions of the BS particles 
(those surrounding strings) become smaller as in Fig. 4b, 



A. Contributions to the mean square displacement 

It is widely believed that the diffusion of a tagged 
or test particle in supercooled and glassy systems is 
caused by its jump or escape motions from temporal 
cages [H E [UHini HMH]. At T = 0.24, we show that 
the diffusion occurs as activation processes, where only 
the BL particles can diffuse over long distances. 

The mean-square displacement M(t) is written as 

M(t) = (|Ar(i)p) ^ l^((Ar,(to,ti))^) (4.1) 

i 

where Ari{tQ,ti) is the displacement length in Eq.(3.6) 
in time interval with width t — ti~tQ. As is well known, 
M{t) exhibits the ballistic behavior (esc t^) in the very 
early stage {t < 1), the plateau behavior M{t) = Mp 
in an intermediate time range, and the linear growth 
M{t) = 6Dt in the late stage. In M{t), the diffusion 
regime starts after dominance of the jump motions over 
the thermal vibrational motions. Since we take the aver- 
age over all the particles in Eq.(4.1), D in this paper is 
the following mean value, 

D = {l-c)Di+cD2, (4.2) 

where Di and D2 are the diffusion constants of the first 
and second species, respectively, and c is the composition. 

In Fig. 6a, the top curve represents M{t), which ex- 
hibits its typical behaviors mentioned above with Mp = 
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FIG. 6: (Color online) Contributions to the mean square 
displacement for T — 0.24. (a) Al{t) from all the particles 
(top curve), Ms(t) from B particles, Mg{t) from BL parti- 
cles, M^{t) from those with An > 0.8, and M(t) for time- 
smoothed positions on a logarithmic scale. See Eqs.(4.4)- 
(4.7). (b) Top five curves represent contributions divided by 
QDt on a semi-logarithmic scale, tending to unity with in- 
creasing t. The approach is very rapid for M^{t)/6Dt and 
Mg{t) /QDt. Shown also is the contribution divided by QDt 
from BL particles belonging to the second species (bottom 
line), which is very close to -D2/(-Di + -D2) = 0.2 for t > 10. 



0.036. Other three curves represent the following contri- 
butions, 



^'^B{t) = ^^((Ar,)20(6,)), (4.3) 

i 

^'^Bit) = (4-4) 

i 
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FIG. 7: (Color onhne) (a) M|(t) vs (^>(t) for various T, 
where both grow linearly in time in the early stage, (b) Ratio 
ifp = Mg{t)/(j)^{t) VS r, where ^jp is a characteristic jump 
length of a particle. 



where Ar^, Bi, and Bf are abbreviations of Ari(io,^i), 
Bi{to,ti), and ;Bf (to,ti), respectively, and = 0.8. 
Here, owing to the step function 6{---), Msit) arises 
from the B particles, Mg {t) from the BL particles, and 
M^{t) from the particles with Ar^ > 0.8. In addition, 
defining the smoothed position 



nit) 



1 

At 



dt'r,{t') 



(4.6) 



t-At 



with At = 100, we calculate the time-smoothed mean 
square displacement, 

M{t)^^J2(\^.(-to+t)-rdto)\'). (4.7) 



In Fig. 6, we can see the sequence MB{t) > M{t) > 
M^{t) > Mg{t). To see the approach to the diffusion 
b_ehavior, in Fig.6b, we display M{t)/6Dt, MB{t)/6Dt, 
M{t)/6Dt, M>{t)/6Dt, and M>{t)/6Dt with D = 
0.596 X 10^^. The latter four quantities approach unity 
more rapidly than M(t) /QDt. In particular, the jump 
contributions M^{t) and Msit)^ become very close to 
QDt from the very early stage t > 50 ~ 5 x IO^'^Tq. 

To be precise, we should calculate Di and D2 sepa- 
rately. To this end, we also consider the contribution 
from the BL particles belonging to the second species, 



M>2(t) = 



N ^ 

ie2 



{{An)'9{B>)). 



(4.8) 



In Fig. 6b, the ratio M^^it) /(^Dt is plotted as a flat line at 
the bottom. It is equal to 0.21 for t > 10 lO^'^r^. This 
value should be equal to CD2/D = D2/{Di + D2), so we 
have the diffusion constant ratio D2/D1 = 0.27 in this 
case. In accord with this ratio, the composition of the 
large particles among the BL particles is c^{t) = 0.25 
for T = 0.24 in Table 1. 

Recall that the fraction of the BL particles 0^(t) 
grows linearly in time in the early stage in Fig. 3, leading 



to the characteristic time r^t in Eq.(3.12). We may define 
a characteristic jump length £jp from the following ratio, 

4 = M>(t)/^>(t). (4.9) 

Now D is expressed as 

D = e%/6T,,. (4.10) 

In Fig. 7a, we plot (j)]^ (t) vs Mg{t) to show their linear 
relationship. In Fig. 7b, we plot £?p vs T, where ^jp is 
close to the particle diameter ai , only weakly depending 
on T. As a natural result, its T-dependence is similar to 
that of the average displacement length (r) bl hi Fig. 4b. 
In our calculation of (f>^ (t) and Mg (t) , the lower bound 
of the jump motions £rn is set equal to 0.8 for any T (see 
Subsec.IVD for comments on this aspect). 

Previously [3 US], M{t) was divided into two contri- 
butions from appropriately defined mobile and immobile 
particles, where the former turned out to dominate over 
the latter in the diffusion regime, in accord with the re- 
sults in this section. 



B. Violation of the Stokes-Einstein relation 

In our model particle system, the violation of the 
Stokes-Einstein relation is ascribed to the steep growth of 
Ta at low T, which has been discussed in Sec.IIIC. How- 
ever, we might need to include other elements (such as 
the nonsphericity of the particles) to explain the observed 
diffusion behavior in real molecular systems [3 [11] [12] . 

In Fig.Sa, we plot D = {Di + D2)/2, D^, and D2 vs T. 
For T > 0.28 we calculated them from the usual mean 
square displacements. Here, D2/D1 = 0.27 for T = 0.24 
and 0.71 for T — 0.40. Thus, D2 decreases more rapidly 
than Di with lowering T, which is consistent with the 
T dependence of the composition of the BL particles in 
Table 1. In Fig. 8b, we plot the products Dra, DiTq, 
D2Ta, , and Dti, vs T. The behavior of -DiTq is nearly 
the same as in the previous simulation [TB]. However, 
DTb depends on T much more weakly than Dto^ where 
Tf, is the life time of bonds in Eq.(3.4). From Eq.(4.10) 
and Fig. 7b, we notice that the relation Tgt = Tb/S well 
holds in the temperature range studied. These results 
indicate that the diffusion is governed by the activated 
dynamics, while the stress relaxation is more sensitive to 
the thermal fluctuations in fragile systems. 

From Fig. 8b, the violation of the Stokes-Einstein rela- 
tion is weaker for a large particle than for a small particle. 
The participation of larger particles in the jump motions 
should become increasingly infrequent with lowering T, 
as suggested by the previous experiments [TTl [T^] . 

C. Van Hove self-correlation function 

The time-correlation function Fs{q,t) in Eq.(2.4) has 
been used so far, though it has been calculated only for 
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0.25 0.3 0.35 0.4 0.25 0.3 0.35 0.4 



T T 

FIG. 8: (Color online) (a) Diffusion constants Di for small 
particles, D2 for large particles, and D = {Di + D2)/2 as func- 
tions of T. (b) D\Ta, D2Ta, and Dt^, and Drb as functions 
of T. The curves of DiTc, and D2Ta indicate the violation 
of the Stokes-Einstein relation from rj/T oc Tc, whereas Drb 
only weakly depends on T. 

q = 27r. Its inverse Fourier transformation Gs(r, t), called 
the van Hove self-correlation function, exhibits thermal 
activation behavior [14i as a function of r = \r\ and t as 
T is lowered. 

From Eq.(2.4) we obtain 

Gs{r,t) = ^Y.(^iAr,{t,,t,+t)-r)), (4.11) 

j 

which is normalized as An drr^Gs{r,t) = 1. It repre- 
sents the probability of finding a particle with displace- 
ment r after time t. In our definition, the average of 
{S{Arj{to,to + t) — r)) is taken over all the particles, so 
it is the mean of the van Hove self-correlation functions 
for the two species as 

G,(r,i) = {l-c)Gsi{r,t) + cG,2{r,t). (4.12) 

Here, Gsaifji) {a — 1,2) are the averages over the par- 
ticles of the species a with N^^J^j Eq.(4.11) being 
replaced by X^jea- The mean square displacement 
M{t) in Eq.(4.1) is expressed as 

/"OO 

M(t)=47r/ dr r'^Gs{r,t). (4.13) 
Jo 

The AI^{t) in Eq.(4.5) follows if the lower bound of the 
above integral is increased to 0.8. 

In Fig. 9, we display 47rr^Gs(r, t) vs r at various t. The 
area below its curve is equal to M{t). For T = 0.24 in 
(a), it exhibits a deep minimum at r = ri„ ^ 0.7. In 
the interior r < it is nearly stationary, while in the 
exterior r > r^ii grows with two or three maxima arising 
from multiple jumps. From Fig. 6, the area of the outer 
region is nearly equal to &Dt for each curve. For T = 0.28 
in (b), in contrast, a shoulder first appears to grow into 
a peak in the exterior r > r,„. The minimum at r = 
increases considerably in time. 
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FIG. 9: (Color online) 47rr*Gs(r, t) vs r at various t. (a) 
For T = 0.24, there is a deep minimum at r = rm ~ 0.7, 
which separates a nearly stationary part for r < rm and a 
growing part with multiple maxima for r > r^- (b) For T = 
0.28, a minimum at r = ~ 0.7 increases in time with one 
maximum for r > r^i- A secondary peak in the outer region 
r > Tui starts to appear as a shoulder. 



From Eq.(4.12), Gs{r,t) is equal to the mean 
Gsi{r,t)/2 + Gs2{r,t)/2 for c = 1/2. We are also in- 
terested in the difference of the motions of the small and 
large particles. Therefore, in Fig. 10, we plot the three 
curves of 'iTrr'^Gs{r,t), 2T:r'^Gsiir,t), and 2Trr^Gs2ir,t) 
at t = 10" for T = 0.24. AU these curves exhibit a 
deep minimum at the common length r = rm- Because 
D2/D1 = 0.27 here, the area of 2Trr'^Gsi{r,t) in the exte- 
rior r > Tin is four times larger than that of 2iTr^Gs2{r, t). 
A second peak corresponding to two jump motions can 
be seen for the small particle but not for the the large 
particle, although the large particles also undergo multi- 
ple jumps on long timescales. 

In Fig.ll, we furthermore divide Gs{r,t) into the con- 
tributions from Bi{to,to + t) = k with k = 0, 1, • • • for 




0.05 



AnrlG, {r,t) 
-27tf^Gs2(r,f) 
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f =10000 



FIG. 10: (Color online) 47rr*Gs(r, t) is decomposed into the 
small particle part 27rr*Gsi(r, t) and the large particle part 
2Trr'^Gs2ir,t) from Eq.(4.12) at t = 10^ for T = 0.24. The 
minimum at r = rm is common for the small and large parti- 
cles. 



three times. The particles with fc > 4 have undergone 
multiple jumps. In (a) at t = 500, those with k < A have 
undergone a single jump, while those with k > 4 two 
jumps. In (b) at t — 10'*, the fraction of the particles 
with fc > 4 is increased. In (c) at t = 5 x lO", the frac- 
tion of three jumps becomes noticeable yielding a small 
third peak. We can thus see how multiple stringlike jump 
motions give rise to these contributions. 

Previously, for the Lennard-Jones potential, Sastry et 
al.fT^ and Wahnstrom [T^ found secondary peaks in 
4Trr'^Gs{r,t) at low T. For hard-sphere binary mixtures 
[TBVllSj . the same behavior was noticed for sufficiently 
large volume fractions of the particles. In particular, Re- 
ichman et a/.|16| found appearance of two peaks in the 
outer region r > indicating discrete step motions of 
hard spheres. 

The behaviors of Gs{r,t) at T = 0.24 can be expected 
generally for glassy particle systems at low T. A parti- 
cle in such systems is analogous to a Brownian particle 
slowly escaping from a cage represented by a potential 
U{r). If the barrier height Um = U{r^) at r = r,„ is much 
higher than fc^T, the probability of staying around this 
barrier is very small (cx exp[— J/m/fesT]) and the proba- 
bility of escaping to the outer region is proportional to t. 
The same behaviors have indeed been found for Gs {r, t) 
aX T — 0.24 in Fig. 9a, where the minimum length is 
analogous to the critical radius in metastable systems. 



D. Dependence on the cut-off length i'm 

In this paper, the lower bound l-^a of the jump mo- 
tions in Eqs.(3.7) and (3.8) has been set equal to 0.8. As 
remarked below Eq.(3.16), this length is between the av- 
erage displacement lengths (Ar)sL and {Ar)BL for the 
BL and BS particles in Eqs.(3.15) and (3.16), where 
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FIG. 11: (Color online) ATvr^G,{r,t) vs r at (a) t = 500, (b) 
10000, and (c) 50000 for T = 0.24. The contributions from 
the particles with Si (to, to + t) = k are shown, which have k 
broken bonds in time interval [to ,to + t]. 



{Ar)BL > 1 and (Ar)BL < 0.5 in Fig.4b. We also note 
that this cut-off length is close to the minimum distance 
0.7 in 47rr^Gs(r, i) in Fig. 9a. In Fig. 12, by setting 
£ni — 0, 7, 0.8, and 0.9, we plot the three correspond- 
ing curves of Mg{t)/6Dt. They surely tend to unity for 
t > 10'^, so Mg{t) is insensitive to a small change of 

At T higher than 0.3, however, the cage life time be- 
comes shorter. As a result, <p%{t) in Eq.(3.10), Tgt in 
Eq.(3.12), and £jp in Eq.(4.9) significantly depend on i^, 
although we have set — 0.8 for any T. 



V. SUMMARY AND REMARKS 

In this paper, simulations have been performed on 
a high-density binary mixture at low T, where the 
particles interact via the soft-core potential (2.1) with 



1 


1 1 

-^AAtf)>0.7 
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7=0.24 
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FIG. 12; (Color online) M>{t)/6Dt from the BL particles 
with r > £m ^ 0, 7, 0.8, and 0.9 for T = 0.24, where = 
0.8 in the other figures in this paper. They are close to the 
minimum distance rm ~ 0.7 of 47rr*G(r, s) in Fig. 9a. These 
three curves converge to unity for t > 1000. 



the size ratio fT2/o'i = 1.2. The total particle number N 
is lO''^ and the density n is 0.8crJ-^ Our main results are 
summarized as follows. 

(i) In Sec. Ill, using the broken bond number Bi(tQ,ti) 
in Eq.(3.5), we have visualized the particles composing 
strings and those surrounding them in Fig.l. The 
former have undergone large displacements (Ar; > 0.8) 
with multiple broken bonds (Bi > 1), while the latter 
small displacements {Ar.i < 0.8) with a single broken 
bond {Bi = 1). The number of the latter is several 
times larger than that of the former because of the 
large coordination numbers in 3D, which leads to the 
difference of the fraction of the former (j>^ (t) and that of 
the latter (j)^{t) = (j)b{t) — (j)^{t) as in Fig. 3. These two 
kinds of motions become increasingly distinct as T is 
lowered as in Figs.l, 2, and 4 and in Table 1. The latter 
contribution to Fs{q,t) in Eq.(2.4) becomes long-lived at 
low T as illustrated in Fig. 5. This leads to the steeper 
rise of than the other characteristic times related to 
the bond breakage at low T in Fig.Sd. 

(ii) In Sec. IV, we have found that the contribution to 
the mean square displacement M(t) from the particles 
with large displacements behaves as 6Dt soon after 
the ballistic regime. The origin of the violation of the 
Stokes-Einstein relation has been ascribed to the steep 
growth of Ta (oc 77) at low T for our system. 

We further give some remarks in the following. 

(1) In our fragile glass-former, the non-Arrhenius behav- 
ior in the Angell plot (the steep growth in the curve of 
log?7 vs 1/T) [2] and the violation of the Stokes-Einstein 
relation [Hi [12] are closely related under the condition 
Ta ^ v/T- Note that we have not calculated 77 but as- 
sumed 77 oc Ttq from Eq.(2.6). In future, we should cal- 
culate 77 together with D for lower T. 

(2) In a 50:50 binary mixture with the soft-core potential, 
we have examined the particle motions around strings. 
We should examine whether our results remain valid or 
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need to be modified for other compositions and for other 
particle interactions [ElESj. We should note that a 20:80 
mixture with a Lennard- Jones potential was used to de- 
tect strings in the original work 'I]. Stringlike collective 
motions were also observed in polycrystalline systems nu- 
merically [24l [25] and experimentally [26j . 
(3) The diffusion behavior of the mean square displace- 
ment M{t) can be obtained soon after the ballistic regime 
if the contribution from large displacements are picked 
up. This is a natural result in the activated dynamics. 
For lower T, the diffusion constant D is obtainable in this 
manner, where its aging behavior will be of interest. 
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